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ABSTRACT 

Context. Theoretical arguments indicate that close-in terrestial exoplanets may have weak magnetic fields, especially in the case of 
planets more massive than Earth (super-Earths). Planetary magnetic fields, however, constitute one of the shielding layers that protect 
the planet against cosmic-ray particles. In particular, a weak magnetic field results in a high flux of Galactic cosmic rays that extends 
to the top of the planetary atmosphere. 

Aims. We wish to quantify the flux of Galactic cosmic rays to an exoplanetary atmosphere as a function of the particle energy and of 
the planetary magnetic moment. 

Methods. We numerically analyzed the propagation of Galactic cosmic-ray particles through planetary magnetospheres. We evaluated 
the efficiency of magnetospheric shielding as a function of the particle energy (in the range 16 MeV < E < 524 GeV) and as a 
function of the planetary magnetic field strength (in the range 0 < M < 10 M^). Combined with the flux outside the planetary 

magnetosphere, this gives the cosmic-ray energy spectrum at the top of the planetary atmosphere as a function of the planetary 
magnetic moment. 

Results. We find that the particle flux to the planetary atmosphere can be increased by more than three orders of magnitude in the 
absence of a protecting magnetic field. For a weakly magnetized planet (M = 0.05 only particles with energies below 512 MeV 
are at least partially shielded. For a planet with a magnetic moment similar to that of Earth, this limit increases to to 32 GeV, whereas 
for a strongly magnetized planet (5Vf = 10.051^0), partial shielding extends up to 200 GeV. Over the parameter range we studied, 
strong shielding does not occur for weakly magnetized planets. For a planet with a magnetic moment similar to that of Earth, particles 
with energies below 512 MeV are strongly shielded, and for strongly magnetized planets, this limit increases to 10 GeV. 

Conclusions. We find that magnetic shielding strongly controls the number of cosmic-ray particles reaching the planetary atmosphere. 
The implications of this increased particle flux are discussed in a companion article. 
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1. Introduction 

1.1. Super-Earths 

Since the first discovery of an extrasolar planet around a Sun-like 
star ( [Mayor and Quelo^ |1995| l, the number of known exoplan¬ 
ets has been growing steadily. Advances in instrumentation and 
analysis have allowed the detection of smaller planets almost 
from year to year. 

The detection of planets with M < lOM®, (i.e., with a mass 
smaller than ten terrestrial masses, which is usually considered 
as the upper limit for super-Earths) became possible about ten 
years after the discovery of the first exoplanets around Sun-like 
stars. The first super-Earth detected by radial velocity measure¬ 
ments is GJ 876d, a planet with ~7.5M® orbiting an M-dwarf 
star (Rivera et al. |2005) l. Other discoveries of planets with simi¬ 
lar masses followed quickly, using a variety of detection tech¬ 
niques. Eor example, OGLE-2005-BLG-390Lb, a planet with 


5.5M(^ in an orbit of 5 AU, was found by microlensing (Beaulieu 
jet al.|[2006) l, and CoRoT-7b, a planet with a radius of 1.7 R® and 
a mass of ~7 — 8M® was detected by transit observations (Leger 
et al. 2009 Queloz et al.j 2009 [ Leger et al. 2011| l. 


The detection limit for planets around main-sequence stars 
has decreased considerably over the last few years: Transit ob¬ 
servations have found planets with radii slightly smaller than 
Earth’s (R — 0.87R®, Eressin et al. , 201^. ap proxim ately half of 


Earth’s radius (R = 0.57R®, Muirhead et al. 


i proxm 

( 2012 |[' 


and even a 


planet with a radius smaller than that of Mercury (R = 0.30R®, 
Barclay et al. [2013 i. At the same time, radial velocity obser¬ 
vations have detected planets with masses similar to that of 
Earth (Dumusque et al. 2012jl. In total, 123 planets with masses 
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evo' 


Dressing and Charbonneau (2013i label this star as potentially 


ved, which would change the derived planetary radius. 
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< lOM® are known to date (November 2014), and 5 planets have 
a mass < IM® 

1.2. M star planets 

One of the many fascinating questions in the held of exoplanet 
studies is the search for habitable worlds. Because of their rela¬ 
tively small mass, low luminosity, long lifetime, and large abun¬ 
dance in the Galaxy, M dwarfs are sometimes suggested as 
prime targets for searches for terrestrial habitable planets (see, 
e.g., [Tarter et al.| |2007| jScalo et ahj |2007| l. For this reason, it 
is interesting to examine super-Earths orbiting M-dwarf stars. 
Currently, 95 of the total of 1850 known exoplanets are located 
around M-dwarf stars (M* < 0.5Mq), and 71 of these are super- 
Earths. While this number is currently still limited, the detection 
of cold debris disks around M-dwarf stars seems to indicate that 
planets may be as frequent for M stars as they are for E, G, and K 
stars (Lestrade et al. [2006 1 . This expectation was conhrmed by 
recent estimations based on the Kepler Input Catalog ( [Dressing 
[and Charbonneau[|2013| l. These estimations indicate that the oc¬ 
currence rate of planets with a radius Q.5R^ < 4/?® orbiting 

an M dwarf in less than 50 days is O.O^qq^ planets per star, and 
the number of planets per star increases with decreasing plan¬ 
etary mass. Thus, the number of currently known super-Earths 
orbiting M dwarfs can be expected to represent only a tiny frac¬ 
tion of the total number. Qualitatively similar results have been 
obtained by Howard et al. ( [2012[ based on KEPLER data) and 
by Bonfils et al. (2013 based on HARPS data). 

A necessary, but not sufficient condition for a planet to be 
considered habitable is the existence of liquid water on its sur¬ 
face (e.g., [Kasting^et^ 1993 Selsis et al. 2007| l. Eor M-dwarf 
stars, the liquid water habitable zone (HZ, the range of distances 
around a star where conditions allow the presence of liquid wa¬ 
ter on the planetary surface) is much closer to the star than for 
Sun-like G-dwarf stars. Its precise location depends on the stel¬ 
lar irradiance (which is a function of the stellar mass), but the 
outer edge of this region typically is <0.3 AU. A first Earth¬ 
sized planet has recently been detected in the circumstellar hab¬ 
itable zone of an M-dwarf star ( [Quintana et al.[ |2014| l, and the 
number of as yet undetected planets at such a location is ex¬ 
pected to be high: Dressing and Charbonneau (20131 estimated 


a mean number of Earth-sized planets in the habitable zone of 
planets per cool star. Similarly, 


M dwarfs of 0.15+[{;^^ 


et al. 


Bonfils 


(j2013|l expect 0.41 ('(q ^3 planets per M-dwarf star. The dif- 


terence between the two numbers is, among other factors, due 
to the different definition of a super-Earth, and the conversion of 
a minimum radius to a minimum mass. The results of [Dressing 
[and Charbonneau (2013' have been reanalyzed with modified 
HZ limits in [Kopparapu ([2013[l, who obtained a frequency of 


0.48('(q 24 and'^AS^^^^terrestrial exoplanets per M-dwarf hab¬ 
itable zone for the conservative and optimistic limit of the HZ 
boundaries, respectively. In any case, all studies agree qualita¬ 
tively: Super-Earths in the liquid water habitable zone of M- 
dwarf stars are very likely to be abundant! 

Care should be taken not to confound this “liquid water hab¬ 
itable zone” with a zone where all planets will indeed be habit¬ 
able. A number of additional conditions needs to be fulfilled for 
“true” habitability (see, e.g., Lammer et al. 2009 |2010[ ). The 
close stellar distance and the special characteristics of M dwarfs 
pose additional problems and constraints to habitability, many 


of which were extensively reviewed by [Tarter et'af] ( [2007[ ) and 
Scalp et ah] (2007 1 . Intense stellar flares are common, especially 
for young M-dwarf stars, but the radiation can be shielded by 
the planetary atmospher e ([Heath et al.[ [1999[ l. Both direct UV 
radiation ( Buccino et al.[ 20071 and indirect UV radiation gener¬ 
ated by energetic photons (Smith et al. 2004[ l can have important 
consequences for biogenic processes on M-star planets. 

A number of additional questions are linked to the presence 
and intensity of a planetary magnetic held. There is increasing 
evidence indicating that super-Earth planets are likely to have a 
weak magnetic held, which has a number of interesting conse¬ 
quences. One of these is the enhanced hux of cosmic-ray par¬ 
ticles to the planetary atmosphere, with potential implications 
ranging from atmospheric chemistry to an increase of the radia¬ 
tion dose close to the planetary surface. 

1.3. Cosmic-ray populations 

The two main populations of cosmic rays are stellar cosmic 
rays (generated by the planetary host star) and Galactic cosmic 
rays (generated by sources outside the planetary system). Eor a 
planet, the Galactic cosmic rays hux can be regarded as isotropic 
and approximately constant (although the hux of low-energy 
particles is slightly modulated by the solar activity); it also only 
weakly depends on orbital distance. Outside the Earth’s magne¬ 
tosphere, it has a peak close to 500 MeV, and the Hux decreases 
for both higher and lower particle energies (Eig.[T] solid line). In 
contrast, the Hux of stellar cosmic rays is time dependent and de¬ 
pends on the stellar activity. Eigure[T] shows three different cases 
for stellar cosmic rays (for a planet at 0.2 AU): The average hux 
during activity minimum (dashed line), the average Hux during 
activity maximum (dash-dotted line), and the cosmic-ray Hux 
during a stellar particle event (dotted line) (see Grenfell et al.[ 
2012[ ?, for details concerning the stellar cosmic ray spectra). In 
either case, the hux of stellar cosmic rays is high for low particle 
energies and steeply decreases for higher energies. At low par¬ 
ticle energies, solar cosmic rays dominate over Galactic cosmic 
rays. Eor Earth, the energy where both contributions have simi¬ 
lar Huxes is between ^100-1000 MeV (depending on solar ac¬ 
tivity). Eor higher energies. Galactic cosmic rays are dominant. 
To compare the relative strength of the two cosmic-ray contribu¬ 
tions at other planets, one has to take into account the planetary 
orbital distance (which has a strong inhuence on stellar cosmic 
rays, but only weakly affects Galactic cosmic rays). Eor exam¬ 
ple, for a planet at 0.2 AU, the energies of comparable Huxes are 
400 MeV, 660 MeV, and 2.6 GeV for our three cases of stellar 
activity (see Eig. [T]). The different energy ranges of stellar and 
Galactic cosmic rays also mean that the particles have different 
penetration depths in the atmosphere, where they encounter dif¬ 
ferent atmospheric species in different concentrations. Eor this 
reason, stellar and Galactic cosmics rays are usually treated sep¬ 
arately. 

The effect of stellar cosmic rays on exoplanetary atmo¬ 
spheres has been addressed by [Grenfell et al.[ ( [2012[ ). An updated 
analysis of this case is presented in a companion article (?). 

In the current work, we focus on the effect of Galactic cosmic 
rays on exoplanets. We will, however, return to the comparison 


between stellar and Galactic cosmic rays in Sect. 4.3 


^ Up-to-date numbers can be found at the Extrasolar Planets 
Encyclopaedia: http: //www. exoplanet. eu jSchneider et al. 201 1|. 


1.4. Cosmic rays on super-Earths orbiting M star 

The severity of the effects generated by cosmic rays depends on 
a number of physical parameters. Previous work has studied the 
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E [MeV] 


Fig. 1. Comparison of cosmic-ray spectra for a planet at 0.2 AU. 
Solid line (black): Galactic cosmic rays. Dashed line (yellow): 
Average of stellar cosmic rays during activity minimum. Dash- 
dotted line (red): Average of stellar cosmic rays during activity 
maximum. Dotted line (blue): Stellar cosmic rays during a stellar 
particle event. 


influence of stellar age (which is a defining variable for the stel¬ 
lar wind parameters, see GrieBmeier et aH|2005| l, the influence 
of the orbital distance ( [GrieBmeier et al.r|2009| l, the influence 
of the presence or absence of tidal locking through its influ¬ 
ence on the planetary magnetic field ( [GrieBmeier et al.| |2005| 
2009 jl and the influence of the planetary type (again, through the 


estimated magnetic held, GrieBmeier et al. 2009| l. Using these 
results, the effects of Galactic cosmic rays on the planetary at¬ 
mospheric chemistry were calculated (Grenfell et al. 2007| l, and 
surface radiation doses were evaluated ( jAtri et al.[|2013| ) Here, 
we take a more general approach and systematically study the 
influence of the planetary magnetic held: Instead of applying a 
model for the planetary magnetic moment, we show how mag¬ 
netic protection varies as a function of the planetary magnetic 
dipole moment. 

Including this shifted focus, the main differences with re¬ 
spect to previous work (GrieBmeier et al. 2005 2009) are the 
following: 


- Here, we systematically study the influence of the planetary 
magnetic held. Instead of applying a model for the planetary 
magnetic moment, we show how magnetic protection varies 
as a function of the planetary magnetic dipole moment. 

- We have included low-energy cosmic-ray particles down to 
16 MeV. 

- We have included high-energy cosmic-ray particles up to 524 
GeV. Thus, the total energy span of the study is (16 MeV 
< E < 524 GeV). The limits of the previous studies were 64 
MeV to 8.192 GeV. We have thus decreased the minimum 
energy by a factor of 4 and increased the maximum energy 
by a factor of 64. 

- The calculation of high-energy particles made it necessary 
to multiply the number of particles by a factor of 4 (to 28 
million particles per configuration and energy). 

- We have extended the analysis of the incoming particle pop¬ 
ulation: We discuss the magnetospheric Alter function and 
compare particle energies and rigidities. We also calculate 
the (magnetic-moment dependent) energy above which a sig- 



Fig.2. Effects of Galactic cosmic rays discussed in this work and 
in the companion article (paper II). 


niflcant fraction of the particles reach the atmosphere, and 
evaluate the energy of maximum particle flux. 


This paper is organized as follows (see also Fig.[^: Section 
[^ describes the types of planetary situations we are interested 
in, namely close-in terrestrial planets around M-dwarf stars with 
weak magnetic fields (Sect. HU- In Sect. |2.2| we present the 
planetary parameters used in our calculations. Sectionj^decribes 
the models and numerical tools we use: The stellar wind model is 
presented in Sect. 3.1 and the planetary magnetospheric model 


is described in Sect. 


3.2 The cosmic-ray model is explained in 


Sect. |3.3| The Galactic cosmic-ray fluxes are analyzed in Sect.j^ 
We have also re-evaluated the implications of such high Galactic 
cosmic ray fluxes on planetary atmospheres and surfaces. These 
implications are discussed in a companion article ([GrieBmeier 


et al.[ 2015 hereafter “paper 11”). Section [^ closes with some 


concluding remarks. 


2. Planetary situation 

2.1. Planetary magnetic moment 


Until recently, the magnetic held of extrasolar rocky planets was 
not only unaccessible to observations, but also to theoretical 
studies. This has changed in recent years, with a number of stud¬ 
ies investigating whether and under which conditions a super- 
Earth can host a significant magnetic held. While these studies 
use very different approaches, they all come to a similar conclu¬ 
sion: Magnetic fields on super-Earths around M-dwarf stars are 
likely to be weak and short-lived in the best case, or even nonex¬ 
istent in the worst case. The relevance of such fields and their po¬ 
tential detectability is discussed elsewhere ( [GrieBmelCT 2014| l. 
With this in mind, the question of planetary magnetic shielding 
against Galactic cosmic rays becomes indeed important. 

In previous studies (GrieBmeier et al. 2005[|2009[ l, we have 
attempted to estimate the magnetic moment of an Earth-like ex¬ 
oplanet around a K/M-dwarf star. Based on analytical scaling 
laws, we found the planetary magnetic dipole moment to lie in 
the range O.OltM^ <M< O.lStM^, where is the value of 
Earth’s current magnetic moment. We then adopted the maxi¬ 
mum value of 0.1551^0 to obtain a lower limit for the cosmic-ray 
flux to the atmosphere. 
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However, such simplified quantitative estimates of magnetic 
fields are not very reliable. More complex approaches, how¬ 
ever, yield values which are not only model dependent, but 
also depend on the precise planetary parameters. For this rea¬ 
son, we take a different approach here: Instead of applying a 
model for the planetary magnetic moment, we show how mag¬ 
netic protection varies as a function of the planetary magnetic 
dipole moment. We thus defer the evaluation of the planetary 
magnetic field, and calculate a large number of representative 
cases instead. This allows us to systematically study the influ¬ 
ence of the planetary magnetic field on the flux of Galactic 
cosmic rays to the planet. In this way, we explore the range 
0.0 M < 10.0 for the magnetic moment, and the 

range of 16MeV < E < 524 GeV for the particle energy. 


2.2. Planetary scenarios 

As mentioned above, previous work has focused on the depen¬ 
dence of planetary magnetic shielding on the stellar age (which 
is a defining parameter for the stellar wind parameters, see 
GrieBmeier et al.[ |2005|l, on the orbital distance ([GrieBmel^ 


et al. 


its ini 


2005 


2009| l, on~the presence of absence of tidal locking through 
uence on the planetary magnetic field ([GrieBmeier et al. 


2009 jl and on the planetary type (again, through the esti 


mated magnetic field [GrieBmeier et al. 20091. Here, we take a 
slightly different approach and treat the planetary magnetic field 

as a free parameter, keeping the other parameters fixed. 

As shown in Table [T] the following parameters are kept 
fixed: The stellar mass M* (M* = 0.5Mq), the stellar radius .R* 
(/?* = 0.46/?q), the stellar age f* (4.6 Gyr), and the orbital dis¬ 
tance d (d = 0.2 AU). These parameters determine the stellar 
wind (which thus remains unchanged, too). In addition, we keep 
constant the planetary parameters of mass Mp (Mp = 1 .OM®) and 
radius Rp (Rp = 1.0/?®). The only planetary or stellar parameter 
we vary in the present study is the planetary magnetic dipole 
moment IM (see Table[T]i. 

We note that the precise values of M*, /?* and of d do 
not matter much for the amount of cosmic rays reaching the 
planetary atmosphere. In particular, using M* = 0.45Mq, /?* = 
0.41/?q and d = 0.153 AU (the values used in the atmospheric 
model of paper II) rather than the values given in Table [T] leads 
to identical particle fluxes within the numerical error. The pre¬ 
cise value of f*, however, can have some limited influence in 
some cases. Based on the results of [GrieBmeier et aL[ ( [2009[ l and 
GrieBmeier et al. (20051, the influence of the parameters d, M*, 


/?*, and f* is discussed in more detail in Sect. 4.2 The main goal 


of this study, however, is to analyze how the flux of cosmic rays 
at the top of the atmosphere varies as a function of the plane¬ 
tary magnetic dipole moment IM and particle energy E over the 
ranges given in Table [T] 

For the planetary magnetic moment IM, we allow values be¬ 
tween 0 and 10 times the terrestrial value (see Table [^. The 
minimum value of the magnetic dipole moment in this study is 
0. This corresponds to an unmagnetized planet, for which we 
take the cosmic-ray energy spectrum outside the Earth’s magne¬ 
tosphere I Seo et al.| ( [1994 1 for £ < 8 GeV, and Mori (1997 using 
their “median” case) for £ > 16 GeV]. The maximum value of 
the magnetic dipole moment in this study is ten times the present 
Earth value, which corresponds to an extremely strongly magne¬ 
tized planet. 

Table [T] also shows the energy range we use for cosmic-ray 
protons. The range of particle energies was chosen such that 


all particles relevant for the processes described in GrieBmeier 
[et al.[([2OT5]l are included (the lowest and the highest energy par¬ 


ticles do not contribute significantly). To serve as input for future 
studies (for which the relevant energy range cannot be known), 
we extended the energy range beyond the values required for 
[GrieBmeier et al.[ ( [2015| l. 


Table 1. Parameters for the different planetary configurations. 


parameter 

value 

M* [Mq] 

0.5"^ 

R* [Rq] 

0.46'' 

u 

4.6 Gyr* 

£/[AU] 

0.2^ 

Mp [M®] 

1.0 

Rp 

1.0 

M 

O.OlM® <M< lO.OlM® 

E 

16MeV < £ < 524GeV 


Notes. M*: stellar mass, /?*: stellar radius, f*: stellar age, d: planetary 
orbital distance, Mp: planetary mass, Rp : planetary radius, M : planetary 
magnetic moment, E: particle energy. 

Notes: ': These values have previously been shown to have little influ¬ 
ence on cosmic-ray shielding (see text). This value has previously 
been shown to have some influence on cosmic-ray shielding (see text). 


3. Models 


In this part, we describe the models used throughout this article. 
In Sect. |3.1| we describe the stellar wind model that we used to 


model the planetary magnetosphere (Sect.[T2]i. In Sect. 3.3 


describe the model used to calculate the flux of Galactic cosmic 
rays through this planetary magnetosphere. The results obtained 
with these models are given in Sect.|^ 


3.1. Stellar wind model 


It is expected that at the close orbital distances of M-star hab¬ 
itable zones the stellar wind has not yet reached the quasi- 
asymptotic velocity regime. Because of the low stellar wind ve¬ 
locity, the planetary magnetosphere is less strongly compressed 
at such distances than one would expect from stellar wind mod¬ 
els with a constant velocity. To capture this behavior realisti¬ 
cally and to correctly describe the flux of Galactic cosmic-rays 
into the atmospheres of close-in exoplanets, we require a model 
with a radially dependent stellar wind velocity. In this work, we 
chose a stellar wind model based on the work by Parker ( 1958[ l, 
which adequately describes the stellar wind around slowly ro¬ 
tating stars. As our study is limited to stars of solar age, this 
approximation is justified. The model takes as input the stellar 
mass, radius, and age. It yields a self-consistent solution for the 
stellar wind velocity as a function of orbital distance v{d), the 
stellar wind density as a function of the orbital distance n{d), 
and the stellar wind temperature T, which we used as input pa¬ 
rameters in the following section. The details of the stellar wind 
model are described in[GrieBmeier et al.[([2005|l and|GrieBmeier[ 
[etar] ( [2009l l. 
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3.2. Magnetospheric model 

The magnetosphere was modeled as a cylinder topped by a 
half-sphere (|Voigt| |1981| |1995[ |Griefimeier et al.| |2004| [2005 [ 


Stadelmann et al.[ 2010 1 . A closed magnetosphere was assumed, 
that is, field lines are not allowed to pass through the magne¬ 
topause. With this model, the magnetic field is defined for any 
point inside the magnetosphere as soon as the planetary mag¬ 
netic moment and the size of the magnetosphere are prescribed. 
As described in Sect. the magnetic moment was varied in the 
range Q.QIM^ < IM < 10.0 fW 0 . The only ingredient missing to 
define the magnetospheric magnetic field thus is the size of the 
magnetosphere. 

The size of the magnetosphere is characterized by the mag¬ 
netopause standoff distance Rs, that is, the extent of the magne¬ 
tosphere along the line connecting the star and the planet. can 
be obtained from the pressure equilibrium at the substellar point. 

This pressure balance includes the stellar wind ram pressure, 
the stellar wind thermal pressure of electrons and protons, and 
the planetary magnetic field pressure: 


2 , O , T- 

mnv^ff + 2nkBT = . 


( 1 ) 


Here, Veff is the effective stellar wind velocity relative to the 

planet (veff = \Jv^ + ''orbit)’ ^ stellar wind ve¬ 

locity, density, and temperature (which can be taken from Sect. 
[^, and /o = 1.16 is the form factor of the magnetosphere and 
includes the ma gnetic field caused by the currents flow ing on the 
magnetopause ( Voigt|[l995[ GrieBmeier et al. 2004 1 . M is the 
planetary magnetic dipole moment. 

We note that Eq. Q does not include a contribution from 
the stellar wind magnetic pressure on the left-hand side. In our 
case (e.g., d = 0.2 AU), v and Bimf are approximately parallel to 
the line connecting the star and the planet. Petrinec and Russell 
( |1997| l analyzed the position of the magnetopause for different 
orientations of the upstream interplanetary magnetic field. For 
magnetic fields Bimf parallel to the stellar wind flow v, they found 
that the substellar standoff distance is not influenced by the mag¬ 
netic field. Hence, in our case, the contribution of the interplan¬ 
etary field to the stellar wind pressure is negligible]^ 

From the pressure equilibrium Eq. ([T]) the standoff distance 
Rs is given by 


Rs = 




(mnvlff -f InksT ) 


1/6 


( 2 ) 


The magnetopause standoff distance thus is a function of the 
stellar wind conditions (obtained in Sect. HD and of the plane¬ 
tary magnetic moment IM. Thus, for each value of the magnetic 


moment, a different magnetospheric configuration is calculated 
according to the cylinder-plus-half-sphere model of |Stadelmann 
et al. (|2010[). With this, the magnetic field is determined for all 


points in the magnetosphere. This magnetic field serves as input 
in the following section. 

3.3. Cosmic-ray calculation 

To quantify the protection of extrasolar Earth-like planets against 
Galactic cosmic rays, the motion of Galactic cosmic protons 
through the planetary magnetospheres was investigated numer¬ 
ically. Because no solution in closed form exists, this is only 
possible through the numerical integration of many individual 
trajectories ( [Smart et al.||200()) l. For each particle energy in the 
energy range 16 MeV < E < 524 GeV and for each magne¬ 
tospheric configuration, over 28 million trajectories were cal¬ 
culated, corresponding to protons with different starting posi¬ 
tions and starting velocity directions. The particles are launched 
from the surface of a sphere centered on the planet with a ra¬ 
dius r > Rs, that is, the particles are launched outside the mag¬ 
netosphere (except for those arriving from the tailward direc¬ 


tion). As compared to our previous studies ( GrieBmeier et al. 
[2005 1 12009| ), this corresponds to an increase of the number of 
numerical particles per configuration by a factor of four. This in¬ 
crease has become necessary because the high-energy particles 
required a better statistical basis. As in any cosmic-ray tracing 
algorithm, the computing-intensive part is not the calculation of 
the particle trajectories, but the evaluation of a complicated mag 
netic field for each particle position. For a specific case. Smart 
et al. ( 2000)1 estimated that the magnetic field calculation takes 


90% of the total CPU time. 

Once the particles enter the magnetosphere, their motion is 
influenced by the local value of the planetary magnetic field. 
This magnetic field is calculated from the magnetospheric model 
determined in Sect. 3.2 using the particle’s instantaneous po¬ 


sition. The trajectories were calculated using the numerical 
leapfrog method ( Birdsdall and Langdon) 1985| l. For each en¬ 
ergy and magnetic configuration, we counted all particles that 
reach the atmosphere, described by a spherical shell one hun¬ 
dred kilometres above the planetary surface (i.e., Rn = Rp-y 100 
km). More details on the numerical calculation of the cosmic-ray 
trajectories can be found in Stadelmann et al.| ( 2010| l. 

The cosmic-ray flux reaching the top of the atmosphere is 
described in Sect. In particular, the flux of particles to the 
planetary atmosphere is quantified by the filter function r|(£’, IM) 
and by the energy spectrum I{E, M), which are defined in Sects. 
4.1 and|4.3| respectively. 


^ There are indeed cases where the stellar magnetic pressure may be 
more important than the stellar wind ram pressure: (1) For very small or¬ 
bital distances, the planet may develop an “ahead shock” in the notation 
of |Vidotto et al.|p010[ ; here, however, we are in the case of a “dayside 
shock”, where the stellar wind ram pressure dominates. (2) One can use 
the magnetic pressure to try to evaluate the size of the magnetopause 
at its flanks rather than at the nose jVidotto et al.||20T3l >. For our case 
(orbital distance of 0.2 AU and using the stellar magnetic field model of 
IGriefimeier et al.||2007[ , this approach yields magnetopause distances 
much larger than the standoff distance (i.e., the pressure balance at the 
nose, using the solar wind ram pressure). To model the magnetosphere, 
we thus have to use the substellar standoff distance, i.e., its size at the 
nose rather than at the flanks, and have to use the ram pressure. This 
shows that for each situation, one has to carefully check which terms 
are relevant in the pressure balance. 


4. Results: Cosmic-ray flux 

We used the numerical model described above (Sect. |3.3] l to cal¬ 
culate the flux of Galactic cosmic rays to the planetary atmo¬ 
sphere. Previously, the cosmic-ray flux to the atmosphere of a 
weakly magnetized {M — 0.15 fW®) Earth-like exoplanet orbit¬ 
ing a K/M-type star with = 0.5M(^ at a distance of d = 0.2 
AU has been calculated ( [GrieBmeier et al.( [2005 [ |2009[ l. It has 
been shown that such a planet will be subject to very high 
cosmic-ray fluxes when compared to Earth. For particle energies 
below 200 MeV, the cosmic-ray flux to the exoplanet was found 
to be up to one order of magnitude higher than on Earth, and for 
energies above 2 GeV, magnetospheric shielding is negligible for 
this exoplanet case. As shown by GrieBmeier et al.[ ( 2009) 1, the 
effect responsible for this reduced shielding efficiency is not the 
compression of the planetary magnetosphere by the enhanced 
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Fig. 3. Magnetospheric filter function r|(£’,f3f). Blue line: case 
of a planet with a magnetic moment identical to that of Earth 
(lAf = l.OfW®). 


stellar wind ram pressure at small orbital distances. Instead, the 
enhanced particle flux is a consequence of the assumed weak 
planetary magnetic dipole moment {fM — 0.1551^0, which was 
estimated assuming tidal locking). Here, we study the question 
how the shielding efficiency varies as a function of the planetary 
magnetic moment. 

For this, we first intro duce the magnetospheric filter func¬ 
tion r|(£’,f3f) in Sect. |4.l[ We then analyze different regimes of 
magnetospheric shielding in Sect. |4.2| Finally, in Sect. |4.3[ we 
evaluate the flux of Galactic cosmic rays to the top of the plane¬ 
tary atmosphere. 

4.1. Magnetospheric filter function r| {E, fM) 

In the magnetospheric model, the magnetospheric filter function 
against cosmic rays is determined in the following way: For each 
value of the particle energy E and planetary magnetic dipole mo¬ 
ment M, the number n,shieided(F, of particles reaching the 
planetary atmosphere is registered. This value is compared to the 
number nunshieided(Fi of particles reaching the atmosphere of 
an identical, but unmagnetized planet. The magnetospheric filter 
function r|(£’, fM) is then defined as 


y\{EM)= (3) 

^unshielded i ^ ) 

Figure [s] shows the dependence of r| on £ and fM based on 
the results of our numerical cosmic-ray model (Sect. |3.3| l. The 
case of a planet with a magnetic moment identical to that of 
Earth (M = 1 .QfMfo) is shown in blue. We define three regimes: 
“well shielded” forr| < 0.10, “partially shielded” for 0.10 < r| < 
0.90, and “unshielded” for r| > 0.90. One can make the follow¬ 
ing observations: 

- For a weakly magnetized planet {fM — 0.05 fWe), the shield¬ 
ing is very weak. At all energies examined in this study 
iE >16 MeV), some protons (at least 30%) may penetrate 
the magnetosphere. Up to 512 MeV, partial shielding pre¬ 
vails. Above 512 MeV, virtually all particles may enter. 

- For a planet with fM = 0.25 cosmic rays below 64 MeV 
are shielded from the atmosphere. Above 4 GeV, all particles 
may enter the atmosphere. 

- For a planet with a magnetic moment similar to that of Earth, 
shielding is almost perfect for £ < 512 MeV. For particle 
energies 512 MeV < £ < 32 GeV, the planet is partially 


shielded. Shielding is virtually absent for energies £ > 32 
GeV. As described in Sect. 4.2 this case is also representa¬ 
tive for the situation on present-day Earth. 

For a strongly magnetized planet {fM = IQ.QfM^), magnetic 
shielding is strong for £ < 10 GeV. For 10 GeV < £ < 200 
GeV, cosmic-ray protons are partially shielded. Only for ex- 
tremly high energies (£ > 200 GeV) is the planet unshielded. 


4.2. Magnetospheric shielding regimes 

To further analyze the result of Fig. we display r| (£, fM) in a 
contour plot. Fig. 4(a) For each planetary magnetic moment iM, 
we define the critical energies E\q%{9v[) and Ei)o%{!M) as the 
energies at which r| = 0.10 and r| = 0.90, respectively, shown as 
solid lines in Fig. |4(a) If we look at the three regimes defined in 
the previous section, we find the planet “well shielded” against 
particles with £ < £io%, “partially shielded” for particles with 
£io% < £ < £ 90 %, and “unshielded” for particles with £ > £ 90 %. 
Thus, the curves defined by Eio%{lM) and Ei)o%{lM) separate 
the three regimes. The observations of Sect. |4.1| remain valid, 
of course, but contours allow a more visual display of the three 
regimes for r|(£,tTf) and of the critical energies Eio%{lM) and 
E‘)o%{fU). 

To compare the behavior of different cosmic-ray particles, it 
is sometimes useful to use their rigidity R instead of their kinetic 
energy £. The kinetic energy £ (expressed in MeV) of a particle 
with rest mass niQ can be converted into rigidity R (expressed in 
MV) using IVogt et aLjdlOOTl Eqs. (11)-(13)): 



£ \ ^ / niQC^ \ ^ 

MeV ) VMeV J 


(4) 


In this, the rigidity R= p/qis the (relativistic) momentum of the 
particle divided by its charge. This choice of variable is useful 
for one specific reason: If one looks at the equation of motion of 
a particle in a magnetic field, the only parameter of the particle 
that enters this equation is the rigidity R. Thus, particles of the 
same rigidity behave in the same way. 

Moreover, for a cosmic-ray particle in a planetary magne¬ 
topause, a change in particle rigidity R can, in some cases, have 
exactly the same effect as a change in magnetospheric configu¬ 
ration, characterized by the magnetic moment IM and the stellar 
wind ram pressure psw ( |Vogt et aT] 2007[ ). This is the case espe¬ 
cially for two locations, namely (a) close to the magnetopause, 
and (b) the inner magnetosphere close to the planet. For these 
cases, we compare two configurations, denoted by subscripts 1 
and 2. The first configuration is thus characterized by fW) and 
Psw.i, and the second configuration by !M 2 and psw, 2 - 


(a) Close to the magnetopause, for two configurations 1 and 2, 
particles have the same behavior if their rigidity Ri and £2 
obey the following relation (|Vogt et al.] |2007| Eqs. (14)- 
(16)): 

Ri \Mj 


(b) In contrast, in the inner magnetosphere, for two configura¬ 
tions 1 and 2, two particles have the same trajectory if their 
rigidity £1 and £2 compare as ( |Vogt et ^|2007[ Eq. (19)): 


£1 _ M 

^ “ W2' 


(6) 
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0.1 1.0 10.0 0.1 1.0 10.0 
magnetic moment (w.r.t. Earth) magnetic moment (w.r.t. Earth) 


(a) Magnetospheric filter function as a function of particle energy (b) Magnetospheric filter function as a function of particle rididity 
and planetary magnetic moment. and planetary magnetic moment. 

Fig. 4. Magnetospheric filter functions t\{E,9v[) and r\{R,9v[) as a function of particle energy E (left panel) or rigidity R (right 
panel) and planetary magnetic moment M. Thick lines: r] = 0.9 and r| = 0.1. Region above r] = 0.9: unshielded. Region between 
T) = 0.1 and Tj = 0.9: partially shielded. Region below r| =0.1: fully shielded. Dashed lines: R°^ M (see text). 


This different behavior in the two regions can serve as a sim¬ 
ple test to see whether the particle shielding is dominated by the 
inner magnetosphere or by the magnetopause. For this purpose, 
it is instructive to reproduce Fig. |4(a)| in terms of the particle 
ridigity R instead of its energy E. Figure |4(b) thus shows the 
magnetospheric filter function r|(R, M) as a function of the par¬ 
ticle rigidity R and the planetary magnetic moment M. Similarly 
to before, the solid lines represent Rio% and R<)q%, that is, the par¬ 
ticle rigidity for which r| = 0.10 and r| = 0.90, respectively. As a 
guide to the eye, the dashed lines obey the relation R M. The 
solid lines should follow the dashed lines when Eq. (j^ holds, 
that is, when the particle shielding is determined in the inner 
magnetosphere. Figure |4(b)| shows that, over the full parameter 
space (e.g., for 0.0514^0 < < lO.Of^g), the solid lines (lines 

of constant Tj) approximately follow the dashed lines. This indi¬ 
cates that the magnetic shielding is dominated by the inner mag¬ 
netosphere, and that the magnetopause contributes only weakly. 


This behavior is not surprising: Inspection shows that 
throughout the magnetosphere, the magnetic field value is close 
to that of a pure planetary dipole, which decreases as B 
This is shown in Fig.|^ which shows the magnetic field (relative 
to the field at the planetary surface) as a function of the distance 
from the planetary center between the planetary surface and the 
substellar point of the magnetopause. The thick solid line rep¬ 
resents a pure planetary dipole. The curved lines represent the 
magnetic field for different magnetospheric configurations (the 
curved line with dots is discussed below in more detail). Close to 
the planetary surface, the magnetic field is identical to the dipole 
case. At larger distances, the field value differs slightly, but the 
difference is small. By construction, the magnetic field value at 
the magnetopause (defined as the location at which the pressure 
equilibrium holds, see Sect. |3.2| i is 2.32 times higher than the 
value for a pure planetary dipole. The locus of these points is 
shown as a dotted line in Fig.|^ The value of 2.32 results from 
the magnetospheric form factor of 1.16 (see references in Sect. 
|3.2| i. However, at that location the magnetic field value is al¬ 
ready considerable lower than close to the planetary surface, so 
that the particle deflection mostly occurs in the inner magneto¬ 


sphere. This explains why in Fig. |4(b)| the particle shielding is 
determined in the inner magnetosphere. 

It is interesting to compare this to previous results for situa¬ 
tions where the orbital distance d and the stellar mass M* were 
varied ( [GrieBmeier et al.||2009| l. Four cases were calculated: d = 
0.2 AU forM* = 0.5 Mq, and B = 1.0 AU forM* = I.OMq, each 
for two magnetic moment cases {M = and M = 9v[^). 

Using the same cosmic-ray model as here, it was found that the 
stellar wind ram pressure (which depends on stellar mass and 
planetary orbital distance) does not have a noticeable influence 
on the cosmic-ray energy spectrum. In other words, the situa¬ 
tion for an exoplanet with d — 0.2 AU for M* = 0.5 Mq is vir¬ 
tually identical to that of a planet orbiting the Sun at B = 1.0 
AU. At that time, this came as a surprise: The stellar wind ram 
pressure does control the size of the planetary magnetosphere, 
thus could have been expected to influence the magnetic shield¬ 
ing of the planet. [GrieBmeier et al.| ( 2009) 1 offered two reasons 
for this behavior; First, for a constant magnetic moment, a de¬ 
crease of the size of the magnetosphere is compensated for by 
an increase of the magnetic field created by the magnetopause 
currents. Thus, a smaller magnetospheric obstacle has a stronger 
boundary, while a larger magnetopause generates a weaker con¬ 
tribution to the magnetic field. This compensatory effect would 
keep the cosmic-ray flux constant when the ram pressure is var¬ 
ied. Secondly, GrieBmeier et al. (2009i also concluded that the 
protons mostly feel the planetary dipole field, and the field of 
the magnetopause currents has little influence. In the light of 
Figs. |4(b) and|^ we confirm this second argument; In the con¬ 
figurations discussed by GrieBmeier et al. ( 2009) 1, the magnetic 
shielding is dominated by the inner magnetosphere, and the mag¬ 
netopause region does not contribute much. The presence of a 
magnetopause modifies the magnetic field value close to this 
boundary layer, but the effect at larger distances is negligible. 
As most of the magnetic shielding occurs in the inner magneto¬ 
sphere, the shielding efficiency is not sensitive to the location, 
or even the existence, of the magnetopause. In other words, the 
protons mostly feel the planetary intrinsic magnetic dipole field 
and the field of the magnetopause currents has little influence 
(cf. [Vogt et al.) [20071 Eq. (20)). This effect alone is sufficient to 
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Fig. 5. Magnetic field between the planetary surface and the sub- 
stellar point of the magnetopause as a function of distance from 
the planetary center, relative to the field at the planetary surface. 
Thick solid line: pure planetary dipole. Curved blue line: case of 
a planet with a magnetic moment identical to that of the Earth 
{!M — 1.0 fW®). Curved black lines: magnetic field for different 
magnetospheric configurations. Straight dotted line: 2.32 times 
the value for a pure planetary dipole (see text for details). 
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Fig. 7. Cosmic-ray flux at the top of the atmosphere I{R, M) as 
a function of particle energy and planetary magnetic moment. 
Dash-dotted line; an unmagnetized planet (fM — 0). Blue line; 
a planet with a magnetic moment identical to that of the Earth 
= 1 . 0514 ). 


explain the results for these configurations. The compensatory 
effect of the magnetic field strength at the magnetopause, while 
valid, probably has a negligible influence. 

Similarly, we can compare our results to those of |Griefimeier| 
et al. (20051. In that work, the stellar age f* was varied. This 


modifies the stellar wind parameters, and thus the size of the 
planetary magnetosphere. In the case of a very young star (f* = 
0.7 Gyr), the cosmic-ray flux to the planetary surface was in¬ 
deed modified, and cosmic-ray shielding became slightly more 
efficient ( [GrieBmeier et al.} |2005| Eig. 8). Our Eig. (curved 
line with dots) shows why: Under the action of the very high 
stellar wind densities and velocities that prevail for young stars, 
the magnetopause moves to very close distances (less than one 
planetary radius above the surface). Here, the magnetopause cur¬ 
rents do indeed modify the whole magnetosphere by increasing 
the local magnetic field strength. Even at the surface, the mag¬ 
netic field strength is increased by ~ 10%, and throughout a 
large part of the magnetospheric volume (approximately half), 
the field strength is increased by at least 40%. In this case, we 
have obviously left the regime where particle shielding is deter¬ 
mined by the inner magnetosphere. Because of the small magne¬ 
tosphere, the magnetopause currents have a non-negligible con¬ 
tribution, and cosmic-ray shielding becomes more efficient with 
increasing stellar wind ram pressure. This has implications for 
the anisotropy of low-energy cosmic ray for a planet around a 
young star. The magnetic field is more compact on the side of 
the planet facing the star than on the opposite side. Thus, the 
flux of Galactic cosmic rays is lower on the star-facing side than 
on the opposite side. Eigure shows that this is indeed verified 
in our simulations (0.7 Gyr case), and that this anisotropy does 
not exist for less compressed situations (4.6 Gyr case). It also 
shows that the impact area is larger when the shielding is less 
efficient (0.7 Gyr case vs. 4.6 Gyr case). 


4.3. Cosmic-ray flux energy spectrum at the top of the 
atmosphere I{E, fM) 

To obtain the cosmic-ray energy spectrum I{E, fM) for extraso¬ 
lar planets, we first established the cosmic-ray flux outside the 
magnetosphere. Eor this, we combined the energy spectra from 
Seo et al. (19941 for £ < 8 GeV, and Mori ( 1997| using their 
“median” case) for £ > 16 GeV. This combined dataset repre¬ 


sents our reference energy spectrum /o(£), which is shown as a 
dash-dotted line in Eig.|^ At the same time, this energy spectrum 
applies to a non-magnetized planet, that is, /o(£) = /(£, = 0). 

We implicitly assumed a stellar environment similar to that 
of the solar neighborhood. The cosmic-ray flux could be differ¬ 
ent, for example, in the case of a different ambient interstellar 
medium ( Scherer et al.[|2002[ [Mtiller et al.j|2006[ [Scherer et al.j 


2006[ 2008|l, or in th e case of the proximity of the Galactic spiral 
arms ( Scherer et al. [[2006) 1. The efficiency of astrospheric shield¬ 
ing (i.e., the fraction of the interstellar cosmic-ray flux penetrat¬ 
ing into the astrosphere) also depends on stellar parameters such 
as the stellar wind ( [Scherer et al.[[2002) [Cohen et al.[[2012| l. 

Strictly speaking, our reference energy spectrum Iq (£) was 
obtained for an orbital distance of 1.0 AU. Not all cosmic rays 
entering an astrosphere reach the inner stellar system. This is 
attributed to effects such as diffusion, convection, adiabatic de¬ 


celeration, and gradient and curvature drifts (e.g., Scherer et al. 
[2006[ [Heber et al.[ [2006j l. Eor this reason, the position of the 
planet within the astrosphere (i.e., the orbital distance) can 
change the particle flux. For distances between 0.1 AU and 1.0 
AU [GrieBmeier et 3T] ( [2009[ Sect. 2.5) found that the particle flux 
varies only weakly (~ 30%), so that we can neglect this effect 
for planets in the habitable zones of K/M-dwarf stars. 

To obtain the cosmic-ray energy spectrum at the top of the 
atmosphere for a given magnetosphere configuration, we mul¬ 
tiplied the magnetospheric filter function r|(£,5Vf) (taken from 
Sect. [4.1[ e.g.. Fig. 0 with the cosmic-ray energy spectrum out¬ 
side the magnetosphere /o(£): 


/(£,5Vf)=Tl(£,5Vf)-/o(£). (7) 

Figure|^shows the resulting energy spectra at the top of the plan¬ 
etary atmosphere (which we considered to be 100 km above the 
planetary surface). The energy spectrum at the top of the atmo¬ 
sphere of a planet with a magnetic moment identical to that of the 
Earth (5Vf =1.0 514) is shown in blue. Figure|^also contains the 
previously studied case of a weakly magnetized ( fM = 0.15 514) 
Earth-like exoplanet orbiting a K/M-type star with M* = Q.SMf^ 
at a distance of d — 0.2 AU ( [GrieBmeier et al.j [2005 [ [2009[ l. 
Figure|^shows how this result is generalized for different values 
of magnetic shielding. A strong magnetic field reduces the flux 
of low-energy particles by more than three orders of magnitude, 
whereas the flux of high-energy particles is barely affected by 
magnetic shielding. 
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(a) Cosmic ray impact region (stellar age: 4.6 Gyr). 


(b) Cosmic ray impact region (stellar age: 0.7 Gyr). 


GrieBmeier et al. 


(2005 I. Left panel: Weakly 


Fig. 6. Cosmic-ray impact region for particles with 256 MeV for the cases presented in 
magnetized exoplanet = 0.15 fW 0 ) at 0.2 AU around an M-dwarf star of 4.6 Gyr age. Right panel: Weakly magnetized exoplanet 
(5Vf = 0.15 Mq) at 0.2 AU around an M-dwarf star of 0.7 Gyr age. Gray regions: No particle impact. 
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Fig. 8. Energy with maximum particle flux to the atmosphere as 
a function of planetary magnetic moment. Middle line: Nominal 
value. Upper and lower line: Estimated uncertainty. 


The most abundant particle population that reaches the plan¬ 
etary magnetosphere is determined by the peak of the curves in 
Eig.|^ The position of this peak (i.e., Fmax = E\i=i^^) is plotted 
against planetary magnetic moment in Eig. Eor magnetic mo¬ 
ments between 0 and lOfW^, the energy of peak flux increases 
from ~300 MeV to 8 GeV. 

Einally, we return to the comparison between stellar and 
Galactic cosmic rays. Stellar cosmic rays are extremely abun¬ 
dant at low energies, that is, below 60 MeV (cf Eig. [^I- Eor a 
planet with !M = 0.25 these particles are strongly shielded 
from the atmosphere. Still, even at energies above 60 MeV, stel¬ 
lar cosmic rays can dominate Galactic cosmic rays: As discussed 
in Sect. 0 for a planet at 0.2 AU, stellar cosmic rays domi¬ 
nate Galactic cosmic rays for particle energies below 400 MeV, 
660 MeV, or 2.6 GeV, depending on stellar activity. As shown 
in Eig. 1^ the magnetic moment required to shield particles of 
these energies are 0.851^0, I.IMq and 2 . 8 fAf 0 (we note that 
these limits are different for planets at different orbital loca¬ 
tions). Eor strongly magnetized planets {M > 2 . 8 fAf 0 ), stellar 
cosmic rays can be entirely neglected, and only Galactic cosmic 
rays matter. Eor planets with an intermediate magnetic moment 
(0.851^0 < M < 2.851^0), stellar cosmic rays are relevant only 
occasionally, for instance, during stellar flares, or during peri¬ 
ods of high activity. For planets with smaller magnetic moments 
(< 0.851^0), both the contribution of stellar and Galactic cosmic 


Fig. 9. As Fig. |4(a)| with shielding of cosmic rays at 400 MeV, 
660 MeV, and 2.6 GeV. Contours: Magnetospheric Alter function 
r|(£’,fTf) as a function of particle energy E and planetary mag¬ 
netic moment 9v[. Region below r| = 0.1: fully shielded. Dashed 
lines: Required magnetic moment to shield cosmic rays of 400 
MeV, 660 MeV, or 2.6 GeV. 


rays matter. As stellar cosmic rays can have a non-negligible 
contribution for planets within a wide range of planetary mag¬ 
netic moments, their influence on exoplanetary atmospheres has 
indeed to be studied; this is done in a companion article (?). 

5. Conclusion 

Magnetic fields on most super-Earths around M-dwarf stars are 
likely to be weak and short-lived in the best case, or even non¬ 
existent in the worst case. With this in mind, the question of plan¬ 
etary magnetic shielding against Galactic cosmic rays becomes 
important. Instead of trying to estimate the planetary magnetic 
moment quantitatively, we calculated a large number of repre¬ 
sentative cases and systematically studied the influence of the 
planetary magnetic held on the flux of Galactic cosmic rays 
to the planet. At the highest energies, we found that the flux 
of Galactic cosmic-ray particles is barely affected by magnetic 
shielding. For lower particle energies, however, we found that 
the particle flux to the planetary atmosphere can be increased 
by more than three orders of magnitude in the absence of a 
protecting magnetic field. Between unmagnetized and strongly 
magnetized planets, the maximum energy for partial shielding 
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increases from 512 MeV to 200 GeV. In the atmosphere, un¬ 
shielded energetic particles can destroy atmospheric ozone and 
other biomarker molecules. Implications include a modification 
of the planetary emission and transmission spectrum and an en¬ 
hanced surface UV flux. In addition, secondary muons can reach 
the planetary surface, leading to a high biological dose rate. 
These effects will be discussed in paper II. For a planet with 
M < 0.25 low-energy stellar cosmic rays can have a strong 
effect that may even dominate Galactic cosmic rays. This case is 
analyzed in a companion article (?). 
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